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ABSTRACT 



In the past decade, helioseismology has revolutionized 
our understanding of the interior structure of the Sun. In 
the next decade, asteroseismology will place this knowl- 
edge into context, by providing structural information 
for dozens of pulsating stars across the H-R diagram. 
Solar-like oscillations have already been detected from 
the ground in a few stars, and several current and planned 
satellite missions will soon unleash a flood of stellar pul- 
sation data. Deriving reliable seismological constraints 
from these observations will require a significant im- 
provement to our current analysis methods. We are adapt- 
ing a computational method, based on a parallel genetic 
algorithm, to help interpret forthcoming observations of 
Sun-like stars. This approach was originally developed 
for white dwarfs and ultimately led to several interesting 
tests of fundamental physics, including a key astrophysi- 
cal nuclear reaction rate and the theory of stellar crystal- 
lization. The impact of this method on the analysis of pul- 
sating white dwarfs suggests that seismological modeling 
of Sun-like stars will also benefit from this approach. 



1. OBSERVATIONAL CONTEXT 



Most of what we can learn about stars comes from obser- 
vations of their outermost surface layers. We are left to 
infer the properties of the interior based on our best cur- 
rent understanding of the constitutive physics. The ex- 
ception to this general rule arises from observations of 
pulsating stars, where seismic waves probe deep through 
the interior and bring information to the surface in the 
form of light and radial velocity variations. The most dra- 
matic example is the Sun, where such observations have 
led to the identification of millions of unique pulsation 
modes, each sampling the solar interior in a slightly dif- 
ferent and complementary way. The radial profile of the 
sound speed inferred from these data have led to such 
precise constraints on the standard solar model that the 
observations and theory now agree to better than a few 

S iarts per thousand over 90 p ercent of the solar radius 
Christensen- Dalsgaardll2002l) . 

If we were to move the Sun to the distance of even the 
nearest star, most of the pulsation modes that we now 
know to be present would be rendered undetectable. We 



would lose most of our spatial resolution across the disk 
of the star, and only those modes of the lowest spheri- 
cal degree (l < 3) would produce significant variations 
in the total integrated light or the spectral line profiles. 
This would reduce the number of detectable modes from 
millions to dozens, leading to a corresponding reduction 
in the ability of the observations to constrain the internal 
structure (e.g., see Kieldsen et al., 1999). Even so, such 
data would still allow us to determine the global proper- 
ties of the star and to probe the gross internal composi- 
tion and structure, providing valuable independent tests 
of stellar evolution theory. 

Recent improvements in our ability to make high- 
precision radial velocity measurements from the ground 
have been driven largely by efforts to detect extra-solar 
planets. These advances in technology have simultane- 
ously led to the first unambiguous detections of solar- 
hke o scillations in other stars (see Bedding & Kields^ 
Eool. Scintillation in the Earth's atmosphere severely 
hmits our abihty to detect the corresponding light vari- 
ations due to these pulsations. But the observational re- 
quirements are similar to programs for detecting extra- 
solar planet transits, so we can expect continued rapid 
progress in this area — primarily from space. 



2. COMPUTATIONAL ASTEROSEISMOLOGY 



Several present (WIRE, MOST) and future (COROT, Ke- 
pler) satellite missions will soon yield nearly uninter- 
rupted long-term coverage of many types of pulsating 
stars. We will then face the challenge of determining the 
fundamental properties of these stars from the data, by at- 
tempting to match them with the output of our computer 
models. The traditional approach to this task is to make 
informed guesses for each of the model parameters, and 
then to adjust them iteratively until an adequate match 
is found. This subjective method is particularly trouble- 
some when combined with a local approach to iterative 
improvement of the defining parameters. 

An optimizatio n scher ne based on a genetic alg orithm 
jCharbonneaul 119951 iMetcalfe & CharbonneauT 12003ft 
can avoid the problems inherent in many traditional ap- 
proaches. Using only observations and the constitutive 
physics of the model to restrict the range of possible val- 
ues for each parameter, genetic algorithms provide a rel- 
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atively efficient means of searching globally for the op- 
timal model. They were inspired by Charles Darwin's 
notion of biological evolution through natural selection. 
The basic idea is to solve an optimization problem by 
evolving the global solution, starting with an initial set 
of purely random guesses. The evolution takes place 
within the framework of the model, and the individual pa- 
rameters serve as the genetic building blocks. Selection 
pressure is imposed by some goodness-of-fit measure be- 
tween the model and the observations. 

iMetcalfd ( 1200 Ih developed a fully parallel and dis- 
tributed hardware/software implementation of the popu- 
lar PIKAIA genetic algorithm written by Paul Charbon- 
neau and Barry Knapp. He used this modeling tool in 
the context of white dwarf asteroseismology, leading to a 
number of interesting physical results including a precise 
estimate of the astrophysically important ^^C(a,7)^^0 
nuclear reaction rate (Metcalfe, 2003|), a nd an empirical 
test o f the theory of stellar crystallization (iMetcalfe et al.L 
l2004h . The impact of this method on the analysis of pul- 
sating white dwarfs suggests that seismological model- 
ing of other types of stars could also benefit from this 
approach. 

We propose to extend this powerful new analysis method 
to treat some of the most important problems in astero- 
seismology, exploiting the full potential of the observa- 
tions. After developing an interface between the existing 
parallel code and models of main-sequence stars, the ini- 
tial applications will be to the well-characterized (5-Scuti 
stars and to the analysis of solar-like oscillations, using 
existing observations and synthetic data. We describe the 
plans for these investigations in greater detail below. 



2.1. Interfacing with Main-sequence Models 

The basic idea behind a genetic algorithm is fairly sim- 
ple: it is just an iterative Monte Carlo method that sam- 
ples the model space randomly, but keeps a sort of mem- 
ory of what worked well in the past. It accomplishes this 
through a computational analogy with the idea of biolog- 
ical evolution through natural selection. It starts just like 
a simple Monte Carlo, where we generate N random sets 
of parameters, evaluate the model for each set, and then 
compare them to the observations. The genetic algorithm 
treats each set of parameters as an individual in a pop- 
ulation, and assigns each a fitness based on how well it 
matches the observations. Next, it selects from this pop- 
ulation at random, with the fittest individuals more likely 
to survive. It then encodes the parameters into simple 
strings of numbers, sort of like chromosomes; it pairs 
them up and performs operations that are analogous to 
breeding and mutation, and then decodes the strings back 
into numerical values for the parameters. This produces 
a new population, so we evaluate the model for each case 
again, and continue the whole process until some termi- 
nation criterion is met. Although it may seem a rather 
contorted way of exploring the model space, there is a 
firm theoretical basis known as the schema theorem to ex- 



plain why it actually works in practice ( lGoldberd.ll989ft . 

Automating the code. Although genetic algorithms are 
often more efficient than other comparably global opti- 
mization methods, they are still quite demanding compu- 
tationally. Fortunately, the procedure is inherently par- 
allelizable; we need to calculate many models, and each 
one of them is independent of the others. So the number 
of available processors determines the number of models 
that can be calculated in parallel. Also, there is very lit- 
tle communication overhead; parameter values are sent to 
each processor, and they return either a list of pulsation 
periods or just a goodness-of-fit measure if the computed 
periods have already been compared to the observations. 
The parallel version of the PIKAIA genetic algorithm 
is perfectly general, and will not require any structural 
modifications to accommodate the application to main- 
sequence models. 

Some customization will be required, however, for the 
main-sequence models that we adopt for this project. 
J. Christensen-Dalsgaard's models were developed for 
the analysis of heli oseismic data, and were used to pro - 
duce 'Model S' of IChristensen-Dalsgaard et alJ (Il996h . 
which has been used extensively as a reference model for 
inversions. Using these models for the analysis of pul- 
sations in Sun-like stars will provide a certain degree of 
internal consistency in our understanding of solar-like os- 
cillations. 

Many of the required modifications will be similar to 
what was necessary for the white dwarf code. The main 
challenge will be to develop a mode of operation that will 
allow the input model to evolve to a specified tempera- 
ture (or luminosity) automatically. A small grid of starter 
models with different masses covering the range of inter- 
est can be constructed independently of the optimization. 
Other interesting parameters, like the helium mass frac- 
tion (Y) and the metallicity (Z), can be specified at the 
beginning of the evolution. When the model has evolved 
to the parameter values requested by the genetic algo- 
rithm, the adiabatic pulsation frequencies can be calcu- 
lated and compared to the observed periods. This will 
lead to a goodness-of-fit measure that the genetic algo- 
rithm will attempt to maximize. 

Optimizing the efficiency. The efficiency of genetic- 
algorithm-based optimization is defined as the number of 
model evaluations required to yield the global solution, 
relative to the number of models that would be required 
for enumerative search of the grid at the same sampling 
density. In practice, a genetic algorithm is usually hun- 
dreds or even thousands of times more efficient than a 
complete grid, and its performance is fairly insensitive 
to the few internal parameters that control its operation. 
We will initially set these internal parameters based on 
our experience with white dwarf models, but we will run 
synthetic data through the optimization procedure (a so- 
called 'hare & hound' exercise) to ensure that the input 
parameters are recovered faithfully. We will repeat these 
exercises with variations in the control parameters until 
the efficiency of the algorithm is optimal. 



3 



2.2. Application to (5-Scuti Stars 

Extending the Cepheid and RR Lyrae instability strip 
down to the main sequence, we find a group of A and 
early-F type stars with pulsation periods between roughly 
half an hour and half a day. This class of pulsating stars, 
with several dozen known members, has been given the 
name of the prototype, (5-Scuti. The amplitude of the ob- 
served variation ranges from a few milli-magnitudes to a 
few tenths of a magnitude, and they often exhibit a mix- 
ture of both radial and non-radial modes. 

Parameter sampling density. The interpretation of pul- 
sation data for (5-Scuti stars is currently facing challenges 
quite similar to those faced by white dwarf modelers 
five years ago. The observational requirements for long- 
term photometric monitoring have been satisfied by suc- 
cessful multi-site campaigns on several stars (e.g., see 
Iregeretal., 1998, 1999;'Handler et al., 2000). But ex- 
ploration of the most important physical parameters in 
theoretical models has been limited to very coarse grids, 
making it difficult to establish a unique best-fit model for 
a particular set of observations. The most extensive at- 
tempt at model-fitting to date, in ter ms of the number 
of computed models, was published bv'Pamvatnvk h et all 
(IT99&) . who explored three parameters with a grid of 120 
models. Using this basic grid, they attempted to interpo- 
late the pulsation frequencies and model parameters to 
produce a much finer grid with 40,000 points. In the 
end they found that the initial grid density was inade- 
quate, leading to significant differences between the in- 
terpolated periods and those resulting from a complete 
evolutionary calculation with the same parameters. 

We will use the genetic algorithm to explore ^-Scuti mod- 
els over a broad range of masses {M) and chemical com- 
positions {Y, Z), covering the full extent of the instabil- 
ity strip and allowing roughly 100 possible values for 
each parameter. With various assumptions about convec- 
tive core overshooting (oov) and rotation (vrot) — or in- 
cluding them as adjustable parameters if warranted — the 
genetic algorithm method will effectively sample the pa- 
rameters with 10-100 times greater density than any pre- 
vious model-fitting attempt. It will do so without the ne- 
cessity of calculating the complete grid, because it sam- 
ples primarily those areas of the model space that it ob- 
jectively finds to produce better fits to the observations. 

Mode identification. The problem with grid density has 
been exacerbated by uncertainties in the identification of 
the spherical degree and azimuthal order (Z, m) of the pul- 
sation modes. The pattern of hot and cool regions on 
the surface of a non-radially pulsating star can be de- 
composed into spherical harmonic functions described 
by these two indices, leading to distinguishable patterns 
of radial velocity and light variation. The frequency of 
the variation also depends on the radial overtone (ri), 
which is not directly observable. The excitation mech- 
anism for the pulsations observed in (5-Scuti stars is not 
well understood; only a small fraction of the pulsation 
modes that are theoretically possible appear to be excited 
to detectable amplitudes. This creates some difficulty 



for a unique interpretation of the observed frequencies, 
a problem which has only recently started to be resolved 
using multi-color photometric techniques ( Dupret et afl 
2003). As space-based data start to become available, al- 
lowing the reliable detection of more low amplitude pul- 
sation modes, this issue is likely to be less of a prob- 
lem. For example, recent WIRE observations of 9^ Tauri 
(Poretti et al., 2002) nearly tripled the number of pulsa- 
tion modes detected in this star, relative to earlier multi- 
site observations from the ground. 

We will investigate the feasibility of incorporating mode 
identification directly into the model-fitting procedure, 
whenever it has not been determined independently from 
observations. The idea is to make no a priori assump- 
tions about the {n, I, m) values in the absence of observa- 
tional constraints, and just try to match the frequencies. 
Within the range of frequencies where pulsations are ob- 
served, the number of modes that are theoretically pos- 
sible increases quickly with the value of /; rotation and 
magnetic fields split non-radial modes into 2l+\ compo- 
nents. So the goodness-of-fit measure must be weighted 
to correct for the mode density, greatly enhancing the fit- 
ness of models that match one of the observed frequen- 
cies with a radial mode (because there are few) and only 
slightly enhancing the fitness when an Z=3 mode matches 
the observations (because there are many). 



Convective core overshooting. Despite these difficul- 
ties, (5-Scuti stars are currently the most promising candi- 
dates for asteroseismology near the main-sequence, and 
objective global model-fitting to their pulsation frequen- 
cies promises to yield important insights into their inte- 
rior structure and evolutionary history. For example, be- 
cause some of these stars are slightly evolved, an exciting 
possibility is to use them as a direct test for the presence 
and degree of convective core overshooting, which is ex- 
pected to leave a dramatic imprint on the pulsation spec- 
trum. Templeton et al. (2001) showed that changes to the 
convective core overshooting parameter in their models 
had a fundamentally different effect on the theoretical 
pulsation frequencies than any of the other parameters 

(M,y,z). 



We will initially perform the model-fitting by considering 
only two possible values for the core overshooting param- 
eter: aov=0.0 (no overshooting), and aov=0.2, the value 
inferred by Ribas et aL ( 2000 ) from an analysis of eclips- 
ing binary data. This will reveal the degree to which the 
frequencies can be fit better by a model that includes con- 
vective overshooting, and it may motivate additional fits 
that allow overshooting to be a fully adjustable parame- 
ter. This is analogous to the way we approached the sen- 
sitivity of white dwarf models to core composition: first 
allowing several discrete values to demonstrate the poten- 
tial of the parameter to improve the fits (Metcalfe et aQ, 
2000), and then following up with a full scale exploration 
leading to an estimate that was much more precis e than 
was earlier thought possible dMetcalfe et alil2001h . 
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2.3. Solar-like Oscillations 

Solar-like oscillations have now been detected in two 
main-sequence stars (a Cen A and B), several sub-giants 
(ry Boo, Procyon, (3 Hyi), and several giants Hya, Arc- 
turus, a UMa). The oscillation amplitudes and the fre- 
quency of maximum power in these stars agree reason- 
ably well with our theoretical expectations. The field 
is progressing very rapidly, and there is good reason to 
believe that many new observations will become avail- 
able in the next few years, particularly after the launch 
of COROT and Kepler With this in mind, it will be a 
distinct advantage to have in place the analysis methods 
that can make sense of these data efficiently, leading us 
quickly to a deeper understanding of the solar oscillations 
in the context of similar pulsations in other stars. 

Extending the method. After adapting the genetic al- 
gorithm fitting method to main-sequence models for the 
analysis of (5-Scuti stars, it will be straightforward to ex- 
tend the method to solar-like pulsators. This primarily re- 
quires a redefinition of the range for each of the relevant 
parameters. We have experience making this kind of tran- 
sition with white dwarf models — extending the genetic 
algorithm method from the helium-atmosphere (DB) to 
the hydrogen-atmosphere (DA) white dwarfs by sim- 
ply changing the allowed temperature range and adding 
one adjustable parameter for the hydrogen layer mass. 
Though not trivial, the time-scale for this development 
was short compared to the time that it took to adapt the 
model to interface with the parallel genetic algorithm. 

A simplifying circumstance for the analysis of solar-like 
oscillations, compared to (5-Scuti pulsations, is the rela- 
tive ease of mode identification. The excitation mecha- 
nism for solar-like oscillations is convection near the sur- 
face, creating a broad envelope of power wi th a peak that 
scales with the acoustic cutoff frequency d Brown et al.L 
1199 Ih . Within this envelope a large fraction of the the- 
oretically possible pulsation modes are excited to de- 
tectable amplitudes, leading to readily identifiable pat- 
terns. Without any detailed modeling, these overall pat- 
terns (the so-called large and small separations, Ai^ and 
5v) immediately lead to an estimate of the mean den- 
sity of the star and can indicate the presence of interior 
chemical gradients. But a full analysis must include a de- 
tailed comparison of the individual frequencies with the- 
oretical models. With a sufficiently long time baseline, 
it should also be possible to resolve the rotationally split 
TO-components of the non-radial (/ > 1) modes. The fre- 
quency separation between these components can yield 
information about the internal rotation rate, since each 
mode samples the interior in a slightly different manner. 

Hare & Hound exercises. Before any new observations 
are available, we will perform theoretical investigations 
using synthetic data to document the results expected to 
emerge from the various space missions, and to bench- 
mark and fine tune the optimization method with 'hare & 
hound' exercises to maximize its efficiency. When new 
observations become available that may benefit from the 
global exploration of models made possible by the ge- 



netic algorithm, we will use this powerful tool to extract 
the complete physical insight that the data can provide. 



3. CALL FOR COLLABORATORS 



The proposal outlined above will initially be supported 
for three years by an NSF Fellowship at the High Altitude 
Observatory, to begin in the fall of 2004. Comments, sug- 
gestions, and collaborators are welcome. Please contact 
Travis Metcalfe <travis@hao.ucar.edu> for more infor- 
mation about the current status of the project. 
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